Orthogonal— Phase— Velocity Propagation of 
Electromagnetic Plane Waves 



Tom G. Mackay^ 

School of Mathematics, University of Edinburgh, Edinburgh EH9 3JZ, UK 

Akhlesh Lakhtakia 
C ATM AS — Computational & Theoretical Materials Sciences Group 

Department of Engineering Science and Mechanics 
Pennsylvania State University, University Park, PA 16802-6812, USA 



Abstract 

In an isotropic, homogeneous, nondissipative, dielectric-magnetic medium that is 
simply moving with respect to an inertial reference frame, planewave solutions of the 
Maxwell curl postulates can be such that the phase velocity and the time-averaged 
Poynting vector are mutually orthogonal. Orthogonal-phase-velocity propagation thus 
adds to the conventional positive-phase-velocity propagation and the recently discov- 
ered negative-phase-velocity propagation that is associated with the phenomenon of 
negative refraction. 
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1 Introduction 

Until quite recently, it was commonly held that the time-averaged Poynting vector and the 
phase velocity of a plane wave in any homogeneous and isotropic medium are coparallel 
[1, 2]; furthermore, the angle between the two vectors was commonly thought to be acute 
in homogeneous anisotropic mediums [3]. The successful fabrication of what are called 
negative-phase-velocity (NPV) materials^ — in contrast to their conventional counterparts, 
the positive-phase-velocity (PPV) materials — led to a radical re-think: the angle between 
the two vectors can be obtuse and even 180°. Much has been written on NPV materials 
during the past five years, and we refer the interested reader to three pedagogical reviews 
[4, 5, 6] and two Focus Issues of electronic journals [7, 8]. 

Could there be materials that are neither of the PPV nor of the NPV types? We have 
already provided an affirmative answer to this question: there can be isotropic dielectric- 
magnetic materials wherein the phase velocity is infinite at some frequency [9]. In addition 

^Corresponding Author. Fax: + 44 131 650 6553; e-mail: T.Mackay@ed.ac.uk. 
^Some researchers prefer the infelicitous term left-handed materials. 
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to the infinite-phase-velocity materials, we provide here a second affirmative answer: mate- 
rials in which the phase velocity is orthogonal to the time-averaged Poynting vector. Such 
materials may be called orthogonal-phase-velocity (OPV) materials. 

The OPV materials emerge here as a gift of the special theory of relativity [10, 11, 12]. 
Relative to a co-moving observer, these materials are supposed to be homogeneous, isotropic, 
and with dielectric and/or magnetic properties. Relative to an inertial observer that is not 
co-moving, planewave propagation may occur such that the OPV condition is satisfied. The 
OPV condition, however, will never be satisfied for the co-moving observer. 

The plan of this paper is as follows: Section 2 is devoted to the Minkowski constitutive 
relations, planewave propagation, and the OPV condition. Section 3 contains numerical 
results for an OPV material that is non-dissipative. Conclusions are presented in Section 
4. In the notation adopted. Re {Q} and Im{Q} represent the real and imaginary parts, 
respectively, of a complex- valued Q. The complex conjugate is written as Q*. Vectors 
are identified by bold typeface and 3x3 dyadics are double underlined; v is a unit vector 
co-directional with v; the unit dyadic is /; and r denotes the spatial coordinate vector. 
The permittivity and permeability of free space (i.e., vacuum) are eo and /Iq, respectively; 
Co = (eoA*o)~^^^ is the speed of light in free space; uj is the angular frequency; and kg — uj/cq. 



2 Analysis 

2.1 Minkowski constitutive relations 

As in a predecessor paper [13], two inertial reference frames, S and E'. provide the setting 
for our analysis. The inertial reference frame S' moves at constant velocity v = vv relative 
to the inertial reference frame E. The electromagnetic field phasors in frames E and E' are 
represented as {D, E, H, B} and {D', E', H', B'}, respectively. The two sets of phasors are 
related via the Lorentz transformation as [10] 

(E • v) V + ^ [(I - vv) • E + v X B] , (1) 
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D' 
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P' 



(l-vv) -B 



V X E 



[(I- vv) -H-v X D] , 
(l-vv) + 



(2) 
(3) 
(4) 



where P = v/cq. 

We consider a homogeneous medium which is characterized by the frequency-domain 
dielectric-magnetic constitutive relations 

B' = HqUj-B' 



(5) 
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in the co-moving reference frame S', with as the relative permittivity and Hr as the relative 
permeability. By virtue of (l)-(4), the constitutive relations (5) transform to the Minkowski 
constitutive relations 

_ _ mv X H 

D = eoCr a • E + 



c, 



•0 



^ mv X E 

B = h iio/ir a • H 

Co - 



(6) 



in the non-co-moving reference frame E; here, 



a — q; / + (1 — a) vv , (7) 
1 

The constitutive relations (6) degenerate to (5) as v — > 0. 

In the remainder of this paper, the propagation of plane waves in a medium described 
by the Minkowski constitutive relations (6) is investigated, with particular emphasis on the 
unusual behaviour which occurs when e^^r < 1- For further general background details 
on planewave propagation within the context of special relativity, the reader is referred to 
standard works [10]-[12]. 

2.2 Planewave propagation 

Our focus lies on plane waves with field phasors 

E = Eo exp {ik • r) 
H = Ho exp (ik • r) 



(10) 



The plane waves are taken to be uniform; i.e., the wavevector k = /ck where the unit vector 
k e and 

k — kR + ikr , (11) 

where kji = Re{k} and kj = lm{k}. 

As a result of combining the frequency-domain Maxwell curl postulates with the Minkowski 
constitutive relations and substituting (10) therein, we find that [13] 

k^k^z^^^;^, (12) 



e(/?k.v)' 



3 



wherein 



A = l + (e,/i, -1)5, (14) 



(/3k. v) 
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>1. (15) 

Since the two wavenumbers represented by (12) are not independent, the medium with 
constitutive relations (6) is unirefringent. The choice of sign for the square root term in (12) 

is dictated by the direction of planewave propagation. 

The amphtudes of the electromagnetic field phasors turn out to be [13] 

Eo = Ciei + 6-2 62, (16) 
Ho = e2 - oueoSr C2 ei , (17) 

where the orthogonal eigenvectors ei and 62 are given by 

ei = kxv, (18) 



62 = a ^ 



f k + m — v) X 61 

Co 



(19) 



whereas Ci and C2 are arbitrary constants. 
Exploiting the relation 

62 = a X ei , (20) 

where 

a^k+ ^'-'-^-' v. (21) 
we find the time-averaged Poynting vector 



|e,pexp(-2a.r) ^ . 

^- ^Re< a + (72 a;eoe„a> 

2 ujfiofi* ) 



(22) 



2.3 Nondissipative mediums 



For nondissipative mediums (i.e., £ M and G 1^), the following three scenarios may be 
identified: 

(i) If EriXr > 1 then A > and the expression (12) represents two real- valued wavenumbers. 

(ii) If < e^/^r < 1 then A can be either positive- or negative-valued depending upon the 
magnitude of /3 and the relative orientations of k and v. Hence, either two real-valued 
wavenumbers or a conjugate pair of complex-valued wavenumbers are represented by 
(12). In the instance where v — 0, we have A = e^/Zr > and the wavenumbers are 
accordingly real-valued. 
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(iii) If Er^ir < then A < and (12) provides two complex-valued wavenumbers. In the 
case where f = 0, we have = which does not correspond to a propagating wave. 
We exclude this scenario from further consideration. 



2.4 Positive, negative and orthogonal phase velocity 

Equation (22) for the time-averaged Poynting vector simplifies to 



eip exp ( —2kj'k • r 



|Ci| 



UJIJ,olJ,r 



+ |C2pa;eoe, 



/cflk + ^ 



(23) 



when Cr G M, /^^ £ and /c e C. The phase velocity 

kR 



(24) 



is described as being positive, negative or orthogonal depending upon its orientation relative 
to the time-averaged Poynting vector. That is, positive phase velocity (PPV) is characterized 
by Vp • P > 0, negative phase velocity (NPV) is characterized by Vp • P < 0, and orthogonal 
phase velocity (OPV) is characterized by Vp • P = 0. In hght of (23), the phase velocity for 
nondissipative mediums is positive, negative or orthogonal depending upon the sign of the 
scalar 



W^kR 



kR + ^(3 (ko-kRpk'v'^k-v 



(25) 



We distinguish between the scenarios A > and A < for nondissipative mediums, as 
follows. 



(a) A > 0. Here the wavenumber is real-valued and 

W ^k[k + ^/3 {ko-k(3 cos 6) cose] , 
with 6 being the angle between k and v. 



(26) 



(b) A < 0. Here the wavenumber is complex- valued, with real and imaginary parts given 

by 

-PC cose 

kR = ko 



kj — ko 



cosey 

1-C{(3 cosef 



(27) 



and W is null-valued. 
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3 Numerical illustrations 



Let us now illustrate the theory presented in the preceding section by means of some nu- 
merical results. In Figure 1, the real and imaginary wavenumbers are plotted as functions 
of the real-valued product 6^//^ for P £ {0.3,0.6,0.9} with 9 e {7r/6, 7r/3, 37r/4}, wherein 9 
denotes the angle between k and v. The wavenumber generally has a nonzero imaginary 
part, despite e and fi being real-valued. Furthermore, the positive-valued upper bound on 
the range of e^/Xr values for which ki ^ increases as the relative speed (3 increases. 

The magnitudes of the phase velocities, corresponding to the wavenumbers presented 
in Figure 1, are plotted in Figure 2. By definition, the phase velocity is ahgned with the 
direction of k. 

The positive, negative and orthogonal nature of the phase velocity, as deduced from the 
sign of the scalar W in (25), are mapped in Figure 2 as functions of (3 and 9 for 

erf^r e {0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5, 2, 3, 4, 5, 10} . 

The OPV region occupies a substantial portion of the P9 plane for small values of er/J^r- As 
Sr/J^r increases towards unity, the OPV region decreases and vanishes in the limit erfj^r 1- 
The OPV region is symmetric about the line 9 = 7r/2, whereas the NPV region is restricted 
to those regions of the (39 plane for which 7t/2 < 9 < tt [13] — which shows that NPV arises 
only when the velocity vector v casts a negative projection onto the wavevector k. When 
erUr — 1, the phase velocity is positive for all values of /3 and 9. For e^/Zr > 1, OPV does 
not arise: the phase velocity is either positive or negative, with the NPV region expanding 
as the product er/ir increases. 

4 Concluding remarks 

Consider a homogeneous medium which is an isotropic dielectric-magnetic medium in a co- 
moving reference frame E', characterized by the relative permittivity scalar > and the 
relative permeability scalar /i^ > 0. From this perspective, the medium supports planewave 
propagation with real-valued wavenumbers, and the phase velocity is directed parallel to the 
time-averaged Poynting vector. 

The nature of planewave propagation is different when viewed in an inertial reference 
frame E, which moves uniformly at velocity — v relative to E'. The following two regimes 
are then distinguished: 

(a) If erUr > 1 then all wavenumbers are real-valued, regardless of v. The projection of 
the phase velocity onto the time-averaged Poynting vector can be either positive or 
negative, depending upon the magnitude and direction of v. 

(b) If < trfJ^r < 1 then planewave propagation characterized by both k & M. and k & C 
(with kj 7^ 0) is supported. Whether kj ^ or not depends upon the magnitude of 
V and its orientation relative to k. For plane waves with kj ^ 0, the time-averaged 
Poynting vector is directed orthogonal to the phase velocity. 
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The regime < e^/x^ < 1 lies between either vacuum (e^ = /i^ = 1) or anti-vacuum 
= fir = —1) and nihiUty (e,. = = 0) [14, 15, 16]. Anti-vacuum is the vaunted material 
for the so-called perfect lenses [17]. Many schemes for materials satisfying the condition 
< ErUr < 1 have been formulated [18, 19, 20], although we must note that dissipation then 
is hard to avoid. However, our calculations (not reported here) show that the OPV condition 
can be very nearly satisfied by weakly dissipative materials. 
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Figure 1: The real (left) and imaginary (right) parts of the wavenumbers fc, normalized with 
respect to fco, plotted against e^/^r for /3 G {0.3, 0.6, 0.9}. The values of are: 30° (solid 
curves), 60° (dashed curves) and 135° (broken dashed curves). 
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Figure 2: The phase velocity, normahzed with respect to Co, plotted against e^/ir for P e 
{0.3,0.6,0.9}. The values of ^ are: 30° (solid curves), 60° (dashed curves) and 135° (broken 
dashed curves). The phase velocity is directed along k. 
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Figure 3: The distribution of positive phase velocity (PPV), negative phase velocity (NPV) 
and orthogonal phase velocity (OPV), as functions of /5 G (0,1) and 6 E (0°, 180°) for 
erfXr = 0.2, 0.4, 0.6 , 0.8, 1, 1.2, 1.5, 2, 3, 4, 5, 10. The labels 'NPV and 'PPV shown hold 
for e,. > and fir > 0, but must be reversed for < and /i^ < 0. 





Figure 3: continued 
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